1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.04 1.13
## dh0 2.13 4.31
## dh1 2.85 4.86
## dh2 1.81 3.05
## dl0 1.45 2.10
## dl1 2.50 4.11
## dl2 1.17 1.47
## dl3 1.82 3.04
##
## Multivariate psrf
##
## 2.54
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1160.44074 | 22017.4891 |
| DIC3 | 1123.74866 | 22087.2419 |
| PWAIC | 26.83821 | 187.8045 |
| WAIC | 1137.43899 | 22099.3759 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.03
## LDevsum 1.01 1.02
## dh0 7.61 14.54
## dh1 8.89 19.10
## dh2 3.34 7.49
## dl0 1.17 1.47
## dl1 1.18 1.47
## dl2 1.03 1.07
## dl3 1.14 1.37
##
## Multivariate psrf
##
## 8.35
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1393.35314 | 27391.3106 |
| DIC3 | 1245.95662 | 27461.9357 |
| PWAIC | 90.50368 | 121.9975 |
| WAIC | 1315.16366 | 27465.8546 |